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^ ' Abstract 

00 ' 

' In this paper a new method is derived for constructing electromag- 

I netic surface sources for stationary axisymmetric electrovac spacetimes 

. endowed with non-smooth or even discontinuous Ernst potentials. This 

' can be viewed as a generahzation of some classical potential theory results, 

, since lack of continuity of the potential is related to dipole density and lack 

■ of smoothness, to monopole density. In particular this approach is useful 

for constructing the dipole source for the magnetic field. This formalism 

^H, involves solving a linear elliptic differential equation with boundary condi- 

I ' tions at infinity. As an example, two different models of surface densities 

^r\| for the Kerr-Newman electrovac spacetime are derived. 

PACS: 04.20.Cv, 04,20.Jb 



1 Introduction 

One of the major challenges of general relativity is the description of a com- 
pact rotating material gravitational source and the vacuum surrounding it. Its 
astrophysical interest is clear, since it would be useful for modelling relativistic 
rotating stars and galaxies. Nevertheless, there are no known exact solutions 
for this problem so far. If we allow the vacuum spacetime to have isometries 
such as stationarity and axial symmetry, the number of solutions of the Ein- 
stein equations that can be obtained by means of generation techniques (cfr. 
for instance [J) is huge. On the contrary there is only a limited number of 
exact solutions, sharing the same symmetries, that could be regarded as inner 
material sources (perfect fluids with a physically reasonable equation of state, 
for instance) and none of these has been smoothly matched to an asymptotically 
flat vacuum spacetime (cfr. j2j for a recent review on the subject). 

In section 12 an outline of the mathematical problem of matching spacetimes 
is provided before dealing with thin layers, that is, two-dimensional sources 
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that could be interpreted as limits of physical configurations in which one of 
the characteristic lengths of the source can be neglected when compared with 
the others. In order to shed light on the method for constructing dipole dis- 
tributions, a reminder of some formulae of potential theory is given in section 
121 These expressions are generalized to stationary axisymmetric electrovacuum 
spacetimes in section 0] The role of the scalar potential is taken on by the com- 
plex Ernst potential j3| , which was introduced for solving the Einstein-Maxwell 
system with a two-dimensional group of isometries on the underlying manifold. 
The lack of continuity of this potential will give rise to a source for the fields 
consisting of monopole and dipole layers. The results will be applied to the Kerr- 
Newman metric in section |31 and compared with the ones achieved by Israel in 
A brief discussion is provided. 



2 The matching problem 

Our aim is the mathematical description of a relativistic compact object and 
the gravitational field in the vacuum that surrounds it. Therefore it is con- 
venient to consider two Lorentzian manifolds (M_|_, respectively 
corresponding to the outer vacuum and the inner material source, whose matter 
content satisfies the Einstein equations, 

Ricci{g±) ~ ^Rg± ^ 87tT±. (1) 

The stress tensor in (M+jg^) can be either electromagnetic or vacuum 
whereas in (M_ ,g-) it can be the one corresponding to an isentropic (charged) 
perfect fluid. 

As a simplification we shall allow both spacetimes to have isometries and 
restrict ourselves to stationary axisymmetric spacetimes. The Killing fields that 
implement these isometries, ^, 77, respectively the generators of the stationary 
and axial symmetry, must fulfil certain conditions. 



5(^,0 <0 g{v,v)>0 [^,v]=0- (2) 

and 77 must have closed orbits. 

The symmetry axis will be then defined by the set of events where the axial 
Killing field is a null vector, 



A = 5(^,^)-0, (3) 

and in order to avoid conical singularities on it, we have to impose a regularity 
condition, namely, |S] 

g(gradA,gradA) 

4A ^ ' ^' 

on approaching the axis. One can also impose that test particles moving along 
the axis do not experience gravitational or electromagnetic forces that push 
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them away from it. This means [H] that in a chart where the coordinates 
, , , x^} are required to satisfy that ^ = d^o, rj = x^ = ^/g{r],r]) 
and that the x^ = constant and x'^ = constant hypersurfaces are orthogonal to 
each other, the connection coefficients vanish on the axis ioi fi ^ 2 ^ i/. 
Furthermore the same restriction on charged test particles implies that the par- 
tial derivatives of the electromagnetic connection along the x'^ direction also 
vanish on the axis. 

We shall also require that not only the metric but also the electromagnetic 
curvature, F, and its four-dimensional Hodge dual, ^*F, and the 4- velocity, the 
density and the pressure of the fluid have zero Lie-derivatives along the Killing 
fields. 

In order to model compact objects such as stars, it is expected that the 
gravitational field will decay to zero at infinity and therefore it is reasonable to 
require that the (M+, (7+) be asymptotically flat in a convenient set of coordi- 
nates. We shall extend on this subject in the following sections. 

Since we are dealing with rotating objects, nonstaticity has to be achieved. 
Otherwise we would have the Schwarzchild solution. Therefore every timelike 
Killing field must be non-surface-forming, 

A dC' ^ 0, (5) 

where b denotes the isomorphism between the tangent and the cotangent bundle 
induced by the metric structure. 

If both spacetimes are to be matched, then there must be in each of them a 
closed 3-dimensional timelike surface E-t that can be imbedded in both g^) 
and {M-,g-) as a submanifold {^±,i±g±) 

z± : S± M±. (6) 

According to Darmois' junction conditions |7] both spacetimes are suitable 
for matching if the induced metric and the extrinsic curvature on the hypersur- 
faces are continuous, 

{^*9]^0 (7) 



[K] = 0, (8) 

after identifying events on S_|_ and E_ by means of an isometry. The difference 
between the values of a quantity on and E_ has been denoted by a square 
bracket. 

If the electromagnetic field F is non-zero, for the sake of continuity it must 
fulfil the condition 

[i*F{n)] = 0, (9) 
on the matching hypersurface, whose outer unit normal is denoted by n. 
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As Israel has shown whenever the equation ||SJ) is not satisfied, the jump 
in the extrinsic curvature K reveals the presence of an energy- momentum surface 
density, S, on E, 

S^^m-Tr[K]i*g), (10) 

and in this case the inner source is surrounded by a 'crust' of matter. 

In the same fashion, there is a non-zero electromagnetic surface current j on 
E, given by the expression, 

[rF(n)] = -4 7rj, (11) 

if the continuity condition Q on F is not satisfied. 

Thin layers are a special case of non-smooth matching. If we try to match two 
vacuum spacetimes across a common hypersurface so that the induced metric 
on both sides fulfils iQ, then the only possibility of having matter content in 
the global manifold is given by (|10|l . (|ll|l . In this case the 'crust' is all we have. 

We would be interested then in drawing physical information about the 
material source that it is located on the hypersurface, such as the densities of 
the physical quantities. In particular dipole densities will be the main concern 
in the following sections. 



3 Classical dipole densities 

Before describing our approach to calculate dipole surface densities in the rela- 
tivistic situation, it will be useful to recall some results of the classical potential 
theory in order to establish a comparison when dealing with general relativity. 

Let us consider a classical field which can be derived from a scalar potential, 
V , that satisfies the Laplace equation. Assume that the field is generated by a 
source located on a closed surface S (a sheet of monopoles and dipoles). Then 
the density of the monopole layer, trg, and the density, (Tm, of the projection of 
the moment in the z direction (its unit vector will be represented by u^). 



dy_ 

dn 



(12a) 



in 



n ■ Uz [V] — z 



dn 



(12b) 



can be calculated in terms of the discontinuities (denoted by square brackets) of 
the potential and its derivative along the outer unit normal, n, to the surface. 

The latter expression can be easily obtained applying the Green identity to 
the flat space fl, excluding the source. 
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0= / (fx{VAz-zAV)^ / dS (v 
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dz 
dn 



dS V 



dz 
dn 



dV 
dn 



dz 
dn 



dV\ 
dn J 
dV 



dn 



(13) 



taking into account that both the cartesian coordinate z and the potential V 
are solutions of the Laplace equation in and that the boundary dfl consists 
of the sphere at infinity and the surface 5*. 

The integral at infinity can be performed if V has the usual multipolar 
expresion, 



V = - 
r 



Q M cos 6* 



0{r- 



(14) 



and it yields the value AttM. Hence the integrand on S can be interpreted as 
the surface density for the dipole moment M of the source. 

This interpretation arises immediately from potential theory (cfr. for in- 
stance [0]), which provides the expressions for the dipole density of a double 
distribution, 



in 



(15) 



in terms of the discontinuity of the potential, and the density of a monopole 
layer, H12a|l . related to the jump of the normal derivative of the potential on the 
surface. 

Hence the first term in (jl2b|) is the contribution of a dipole layer to the total 
dipole density and the second term is just the moment density in the z direction 
of the distribution of monopoles. 



4 Relativist ic layers 

For the stationary axisymmetric Einstein-Maxwell system we shall follow a sim- 
ilar approach, which can be viewed as a generalization of the formalism intro- 
duced in ^U], 22 for static electrovacs. In that paper magnetic surface-sources 
were constructed for magnetostatic solutions of the Einstein-Maxwell equations. 
The magnetic moment density was calculated from the discontinuities of the 
Ernst magnetic potential involved in the generation of the solutions. How- 
ever that formalism could not deal with asymptotically monopolar electric fields 
and with nonstatic metrics and, therefore, many important exact solutions fell 
out of its scope. 

In order to work with the Ernst formalism for the stationary axisymmetric 
Einstein-Maxwell system, we shall work in charts adapted to the Killing fields 
throughout the paper so that these can be expressed in the form, ^ = 9t, ?7 = 9^, 
and restrict ourselves to metrics which can be rendered into a reducible matrix in 
a holonomic frame, that is, there is a set of coordinates in which the metric can 
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be expressed as g = gi® 92, where gi is the metric in the subspace spanned, at 
each event, by the Kilhng fields and 172 is the metric in the subspace orthogonal 
to the former one. Having this restriction in mind, it is useful |12) to introduce 
an orthonormal frame, {0°, 6^, 9^, 9^}, such that the Killing fields have non-zero 
projection on 0° and 6^ only. 

In this orthonormal frame, we shall consider electric and magnetic forms, E 
and B, with zero projection on the orbits of the Killing fields, 

F =^ E A9" + AO^ ^ *F = -B A9^ + ^E A9^, (16) 

so that the electromagnetic Faraday, F, and Maxwell,'' * F, forms have a simple 
expression after introducing the Hodge duality (*6'^ ~ 9^) in the space orthog- 
onal to the Killing fields jlOj. 

Taking into account Cartan's structure equations, expressing the torsion- 
free connection coefficients as one-forms, a, w, s, b and (cfr. ^HIiCSI for the 
physical meaning of these forms) 

d9° = aA9^ -w A9^ (17a) 
dO^ = {b-a) A9^ - s A0° (17b) 
d9'^ ^-i'A9^ (17c) 

de^ = i'A0'^, (17d) 

the Maxwell vacuum equations, dF = and d * F = 0, can be written in 
a compact fashion after defining a complex one-form f = E + iB, 

df ^~aAf ~isA*f (18a) 

d^ f = {a-b) A^f + iw A f. (18b) 

The orthonormal frame can always be chosen without loss of generality so 
that the one- form s is zero. Bianchi's compatibility conditions for ()17f) . 

db^O (19a) 

da^w A s (19b) 

dw = -{b-2a)Aw (19c) 

ds^{b-2a)As, (19d) 
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can be then formally integrated after introducing some new functions U, A, 

a^dU b^dlnp (20) 



w = p-^e^^-'dA s = 0, (21) 

which allow the solution of Cartan's first structure exterior system of equations 
p7|) and yield the line element in Weyl coordinates, 

ds^ = -e2^(dt - Ad0)2 + e-2^{p2 rf02 ^ g2fe(^^2 ^ ^^2)1^ (22) 

in terms of the functions U, A, k oi p and z. 

Also equation H18a(l can be solved in terms of a complex scalar potential, <i>, 
the Ernst electromagnetic potential 

f ^ E + iB ^ -e^'^ d^, (23) 
that satisfies one of the Ernst equations , 

d*d<^ + {b-2a) A*d<^ = i dAAd<^, (24) 

which can be easily obtained from equations IjlSbp and (|23|) . However, it will 
be convenient to cast it in a different form, 

d(e"^'^p*d$-iAd$) = 0, (25) 

for future purposes. In a coordinate patch, it bears a resemblance with a complex 
Laplace equation on the hypersurfaces of constant time, 

L$ = -^a,, |iVV5 (^e-^^ g'^'^ -'-Ae'^"^ 5,$|-0, (26) 

but including a correction depending on A due to non-staticity, that prevents 
the decoupling of the equations for the real and imaginary parts of $. In this 
equation g is the metric induced by (|22|l on the hypersurfaces t = const., N = 
(— '*(7**)~2 is the lapse function and e is the Levi-Civita tensor on the surfaces 
of constant time, t, and azimuthal angle, (f). For simplicity the whole equation 
has been written as the action of a differential operator, L, on the potential. 

As a matter of fact, this is a consequence only of the Maxwell equations 
in the curved spacetime whose metric is given by (I22II . regardless of whether 
the electromagnetic field is the source of the gravitational field, since we have 
not made use of the Einstein equations. Therefore what follows is valid for 
any potential that satisfies the equation L$ = on the previously described 
geometry. 

As it was stated in the first section, we are interested in compact sources and 
therefore we shall only consider metrics which can be rendered asymptotically 
flat in some coordinates {t, r, 6, 0), 



7 



, 9 / 2m, , , 2Jsin^ ^ , ,x9 

ds^ = -1 )idt+ d(j))^ + 

r r 

+ (1 + — ){dr2 + (r2+cir)(d02^sin2 0d02)} + O(l/r2), (27) 
r 

that allows us to read the total mass, m, and the total angular momentum of 
the source, J, from the Lense-Thirring expansion. In this set of coordinates, 
the electromagnetic potential of the compact source will be required to have an 
asymptotic expansion in terms of the first multipole moments, the total charge 
e and the complex electromagnetic dipole A/, whose real and imaginary parts 
are, respectively, the electric and magnetic dipole moment, 

^ e M cos 9 C2 ^, Q, , 

The constants ci and C2 may arise in some choices of coordinates. 

In the case of a non-smooth or discontinuous electromagnetic Ernst poten- 
tial, $, the equation i"I> — may not hold in the whole manifold, and it may 
have a distributional source located in the region where the discontinuity takes 
place. Since the aim of this paper is devoted to thin layers in stationary ax- 
isymmetric spacetimes, we shall assume that this region is a closed surface S of 
outer unit normal n within the constant time hypersurfaces. We shall denote 
these hypersurfaces as (Va,^) for our calculations. 

In order to mimic the construction done in section it will be necessary to 
obtain a Green identity for the L operator. A most natural candidate can be 
checked to be, 



/ ^{ZL^-^^L+Z)dx^dx^dx^ ^ 

= [ dS N U-^" ( ^^^] +i- (Z *d^(n) + ^ *dZ(n))\ , (29) 
Jan I \ dn dnj p J 

which is a consequence of the divergence theorem. i+ is just the complex 
conjugate of L. We shall apply this identity to the region O = U Vg", the 
disjoint union of the outer and inner regions of V3 referred to the surface S. The 
oriented boundary of is just the surface S whereas that of is formed by 
the sphere at infinity and S. 

So far no condition has been imposed on the function Z. In the classical 
case, it was just the cartesian coordinate z and the Laplace equation held for it. 
Since we no longer have a Laplace operator but L, it seems natural to choose 
Z so that it satisfies L~^Z = in 17. Hence the left hand side in H29|l will be 
zero and we shall have to deal with surface integrals only. In addition to the 
differential equation the function Z will be required to behave near infinity as 
the cartesian coordinate. The set of conditions on Z will be then, 

L+Z = Z= (r + cg) cos 61 + O(r-i), (30) 
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which requires solving an elhptic partial differential equation with boundary 
conditions at infinity. Again C3 is a constant. 

The integral at infinity in H29|) can be calculated from the information given 
by the asymptotic behaviour. We are left with just an integral over the surface 
S, where the source is located, 



= -AttM+ I dSN {e-^^ 



$ ^ _ Z — 

dn dn 



A 



i — [Z * d<P{n) + $ * dZ{n)] 
P 

(31) 

which allows us to express the total electromagnetic moment as an integral over 
the source, 



M = — / dSN 

4tt 



-2U 



^dZ 

dn dn 



-i—[Z* d<^{n) + $ * dZ{n)] 
P 



As it was done in the classical potential theory, we can interpret the discon- 
tinuity of the integrand as the electromagnetic dipole moment surface density 
of the source for the field. 



47r 



■N le 



-2U 



dn dn 



i-\Z * d$(n) + $ * dZin)] \ , (33) 
P 



the real part of ctm being the electric dipole density and its imaginary part, the 
magnetic moment density. 

When A is zero, that is in the static case, this formula is pretty similar to 
the classical one (jl2b|l . corrected by metric factors. Its first term also coincides 
in the static case with the formula introduced in 10 , where the contribution of 
the moment density arising from the layer of monopoles (second term) was not 
taken into account. 

Also the electric charge density on S can be calculated using the divergence 
theorem, since equation H26|l has the form of a total derivative and therefore its 
integral on f2 can be reduced to a surface integral on its boundary, 

0= / L$= / dSN[e-^^'^+i- *d<^{n)\ . (34) 
Jan I dn p J 

The asymptotic expansion of the fields provides the necessary information 
to perform the integral at infinity, 







-47rQ 



dSN { e~2^ 



dn 



i — [*d$(n)] 
P 



(35) 



from which we can read the expression for the charge density, 



Q 



dS aq 









dn 



[*d$(n)]^ I , (36) 
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in terms of the discontinuities of the derivatives of the electromagnetic Ernst 
potential This formula recovers Israel's expression for the electric charge 
density but written in terms of the electromagnetic Ernst potential. 

The function Z can only be considered as a coordinate in the static case. 
Otherwise the equation H3U|) is complex and so is the solution Z . This is very 
similar to what happens with the Ernst potential 0,2 = + *X- Whereas 
its first term is a metric function, the second term is just an auxiliary potential, 
the twist potential x, that is due to nonstaticity. Similarly, just the real part 
of Z can be viewed as the coordinate function which determines the projection 
of the dipole moment that is being calculated. The imaginary part is again an 
auxiliary potential which states the influence of the electric (magnetic) field on 
the magnetic (electric) dipole density. 

As it has already been mentioned, this results are valid for any potential that 
satisfies the equation in the geometry defined by Therefore the same 
considerations can be applied to the gravitational Ernst potential e P] , since it 
fulfils Le = 0. It will be convenient, however, to introduce another potential 77, 



r; - e - 1, (37) 

since the Ernst potential is defined up to a constant, which is usually fixed so 
that the potential tends to one at infinity and this asymptotic finite value of e 
is not very convenient when performing integrals on fJ. If the potential has the 
following asymptotic behaviour, 

2 TO 2 J cos 6* C4 ^, o, 
^ = -— - * + + 0{r-\ (38) 

where C4 is again a constant, then all the previous calculations can be repeated 
just substituting <I> for r\ to yield expressions for the mass and angular momen- 
tum densities, 

+ z- [*rfr/(n)]l (39) 
P J 



= — TV 
Stt 



-2(7 



dn 



on the surface S 



-2U 



^ dr/ dZ 
dn ' dn 



i — [Z * difln) 
P 



f] * dZ(n)] 



(40) 



5 An example: The Kerr-Newman spacetime 

As an application of this formalism, we shall calculate the sources for the phys- 
ical quantities of the Kerr-Newman spacetime ^Sl- The energy-momentum 
tensor of a source located on the r = surface of this spacetime has already 
been obtained by Israel |4j using the formalism of thin layers |3 and afterwards 
by Lopez ^3] using distributions. Wc shall focus our attention on the magnetic 
moment density since it is not calculated in those references. 
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Instead of allowing the radial coordinate r to become negative, we shall 
restrict its range to positive values, as it is done in j^l- This amounts to iden- 
tifying points on the hypersurface r = (both sets of coordinates {t, (p, 0, 6) 
and {t, 0, TT — 0) represent the same event), as if we were working with oblate 
spheroidal coordinates. This identification causes that certain functions of the 
collatitude angle will be discontinuous, such as the cosine, and hence the Ernst 
potentials, 

2 m e , , 

r 1 a cos r — t a cos 

and their derivatives will encounter discontinuities on r = 0, revealing the pres- 
ence of the source. 

The Kerr-Newman metric, in Boyer-Lindquist coordinates. 



, ., , ,,,, (2mr ~ e^)a sin^ 



r 



a? cos^ 9 j-2 _ 2mr + + a? cos^ ( 



+ (1 - ^"^^ - 2mr + + e^) sin^ 0d^^ + 

^ +a? cos2 0' ' 

+ {r^ -2mr^c?co^^0^e^){- ^ + (10"^)} (42) 

if the parameters m, a, e satisfy < -|- e^, induces a line element on the 
surface r 0, i =const. 

dsi = cos^ d0^ + sin^ 6* (a^ - tan^ 6*) dcp^ , (43) 
and unit normal, 



n=^^a., (44) 
a cos 6' 

where m is the total mass, ma, the total angular momentum and e, the total 
electric charge. 

It is clear that the Kerr-Newman metric and Ernst potentials fulfil the pre- 
viously described asymptotic requirements and therefore the surface densities 
can be calculated within this formalism. We need a solution for 13U|I . 

r7 / o \ /, COS0 + i amcos^ ^ 

Z = (r - 2to) cos6' H , (45) 

^ ' r + ia COS0 ' 

before calculating the densities on the surface r = by introducing these ex- 
pressions in (|^ . ijTO . 

From the expressions for the electromagnetic moment and the charge surface 
density, 

(e^ cos^ + e'^ + a'^ cos^ 0) ie , , , , ^ , , 

= — e e [0, 0o) U (00, ^/2), (46a) 

2 TT a-^ COS'' 6^ — tan"^ 0\ 2 
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^Q = -- TT-, T—Y^ 0e [O,0o)U (00,^2), (46b) 

2 TT a cos'^ 0\a'^ — tan \ 2 

where Oq = arctan(|a/e|). The corresponding angular momentum and mass 
densities are obtained by muhiplying them, respectively, by the inverse of the 
gyromagnetic ratio, e/m, since the Ernst potentials are linearly dependent, 



(e^ cos^ 9 + + 0? cos^ 6) im , . , , . , . 

= T^2 2. /flii ^ e [0' U (^0, 7r/2), (46c) 

2 TT cos'^ a-^ — tan 2 



*T1 

[O,0o)U(0o,7r/2). (46d) 



2 TT a cos-^ 9\a'^ — tan^ 6*1 2 

The electromagnetic dipole density, aM^ is imaginary, and hence there is no 
electric dipole density, just the magnetic moment density. 

From the expressions for the differential elements of magnetic moment and 
electric charge. 



(e^ cos^ 6* + + cos^ 0) e sin^ , , 

dM^ag = 1{oM } dS = ^ — '- dQd4>, 47 

2 vr a COS"' 9 

dQ^agdS^- /'"'"^^ de d(t>, (48) 

2 TT COS^ 

we could get by integration on the surface r — the total magnetic moment 
and electric charge of the source. However, instead of obtaining the values 
Mjnag = ea and Q = e, the integration yields divergent results due to the 
singular ring at r = 0, 6* = n/2. More precisely the integrals 

i'2tt /"aiccos ^ / \ 

M = lim / d(h i dOoM ^iea{\ + lim — | (49a) 

Q = Urn y d(P J d9aQ= e (l - lim - j , (49b) 

provide the correct results plus a term that blows up on approaching the 
singular ring. There is an infinite contribution to the physical quantities located 
on the ring. 

It is remarkable that the angular momentum density is integrable in the 
uncharged case, as it was shown in jl5|. From the expression of the induced 
metric (I43f) we learn that there is a change of signature at da. For values of 6 
greater than the metric on the surface is semiriemannian. 

The elements of mass and electric charge are the same that were obtained 
in ^ using the thin layer formalism. The linear relation between the magnetic 
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surface density and the angular momentum surface density confirms the value 
of the gyromagnetic ratio, e/m, for the Kerr-Newman spacetime. 

There is another choice of location for the source of the Kerr-Newman space- 
time, as it is shown in JS], J7]. Instead of considering the surface r = 0, it is 
possible to match the Kerr-Newman manifold to Minkowski spacetime at the 
pseudosphere r ~ vq ~ j^m. The electromagnetic Ernst potential in the inte- 
rior of the pseudosphere is taken to be zero. An advantage of this choice is that 
the source is valid for all values of the parameters m, a, e, since this surface 
is always located in the Boyer-Lindquist chart (the horizon, whenever there is 
one, lies within the surface). Another advantage of this approach lies on the 
fact that the metric at r = tq. 



ds^ = -dt^ + {rl + cos^ B) d0^ + — + (r^ + a^) sin^ 6 dc/)'^, (50) 

V ^0 + / 

is just flat spacetime in oblate spheroidal coordinates, as it happened for the 
uncharged case at r = 0. Therefore the metric tensor is continuous at the 
matching surface and also the induced metric, 

dsl = (^0 + cos^ 6) dO^ + (r^ + a^) sin^ 9 dcj)^, (51) 

on the pseudosphere experiences no change of signature. 

The expressions that we obtain for the electromagnetic dipole and monopole 
densities on the surface r = tq. 



e cos 



9 {2rl {ro-m) + ia cos^ (a^ cos^ 9 + 3r^)} \/rfT^ 



4 TT (rQ + cos^ 



(52a) 



a' cos' 9) \/r^ 



- \ , ^ ^ , (52b) 

47r (rg + cos^ 0) ' 

lead to the correct results respectively for the total magnetic moment and 
electric charge after integration on the pseudosphere. From the fact that gm is 
not imaginary, we learn that there is a non-zero density of electric moment on 
the matching surface, given by the real part of (|52a|l . 



2 e cos 6* rl (ro - m) .Jrl + 
4 7r (r2 + a2 cos^ ef/^ 

although the total electric moment amounts to zero. The magnetic moment 
density is then given by the expression 



e a cos^ 9 (a^ cos^ 9 + Stq) y^r^ + a'^ 
^A^.-, = ; — — TTT^T^ • (54) 



4 7r (r2 
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6 Discussion 



In this paper a new method for constructing dipolc surface sources for stationary 
axisymmctric clcctrovac spacetimes has been introduced. The expressions that 
are obtained show the contributions of dipole and monopole layers to the total 
dipolc density. For this purpose it is necessary to calculate a function Z as 
the solution of a linear elliptic differential equation with boundary conditions 
at infinity. Its real part can be interpreted as the function which determines 
the projection of the dipole moment that it is being calculated whereas the 
imaginary part is an auxiliary potential related to nonstaticity. 

The method has been applied to the Kerr-Newman spacetimes to produce 
the magnetic moment density of two different choices of source, both of them 
on hypersurfaces of constant time. To our knowledge, this quantity had not 
been calculated before. One of them has been located, for parameters m, a, 
e satisfying w? < + e^, on the surface r = 0, that is surrounded by a 
singular ring, whereas the other lies on the pseudosphere ro = e^/2m. The 
resulting densities are integrable in the pseudosphere model but they are not 
so in the first model. Besides the pseudosphere model is valid for all ranges of 
the parameters m, o, e, since the horizons lie within the surface. Therefore the 
second choice seems to be a more suitable source for the Kerr-Newman fields. 

For the future it would be interesting to devise new methods for interpreting 
the lack of continuity or smoothness of the potentials in terms of distributions 
in order to cope with other possible sources, such as struts and rings. 
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and by a FPI Predoctoral Scholarship from Ministerio de Educacion y Ciencia 
(Spain). The author wishes to thank F. J. Chinea and L. M. Gonzalez- Romero 
for valuable discussions and the referees for their suggestions. 
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